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Abstract 

A two-state spin system is specified by a matrix 
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where /3, 7 > 0. Given an input graph G = (V, E), the partition function Z\{G) of a system is defined as 

Za(G)= ]T II A *(u)Mv)- (!) 

cr:V-^{0,l} («,t))£E 

We prove inapproximability results for the partition function in the region specified by the non-uniqueness 
condition from phase transition for the Gibbs measure. More specifically, assuming NP 7^ RP, for any 
fixed p, 7 in the unit square, there is no randomized polynomial-time algorithm that approximates ZpSG) 
for rf-regular graphs G with relative error e = 10~ 4 , if d = Q(A(/?,7)), where A(f3 7 j) > 1/(1 — $7) is 
the uniqueness threshold. Up to a constant factor, this hardness result confirms the conjecture that the 
uniqueness phase transition coincides with the transition from computational tractability to intractability 
for Zx{G). We also show a matching inapproximability result for a region of parameters /?, 7 outside the 
unit square, and all our results generalize to partition functions with an external field. 
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1 Introduction 



Spin systems are well studied in statistical physics and applied probability. We focus on two-state spin sy- 
stems. An instance of a spin system is a graph G = (V, E). A configuration a : V — > {0, 1} assigns to every 
vertex one of two states. The contributions of local interactions between adjacent vertices are quantified by 
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a 2 x 2 matrix with /3,7 > 0. The partition function Z A (G) of a system is defined as 

Za{G) = ^ Y\ A a(u),a(v), (2) 

<t:V-s.{0,1} (u,v)GE 

and we use uj(G, a) = JT^ U V ^ &E A a ^^^ to denote the weight of a. 

For a fixed A, we are interested in the complexity of computing Z A (G), where G is given as an input. 
Many natural combinatorial counting problems can be formulated as two-state spin systems. For example, 
with f3 = and 7=1, Z A (G) is exactly the number of independent sets (or vertex covers) of G. The defi- 
nition of Z A (G) in ^ can be generalized to larger A, and the problem is also known as counting (weighted) 
graph homomorphisms [20^ I16j . On the other hand, the Ising model is the special case where f3 = 7. 

The exact complexity of computing Z A (G) has been completely solved for any fixed symmetric matrix 
A [TTJ [3l [131 E] and even for not necessarily symmetric A [8], HI El E] as part of the dichotomy theorems 
for general counting constraint satisfaction problems. When specialized to two-state spin systems, Z A (G) 
is #P-hard to compute exactly, except for the two restricted settings of f3j = 1 or (3 = 7 = 0, for which 
cases it is polynomial-time computable. Consequently, the study on two-state spin systems has focused on 
the approximation of Z A (G), and this is the subject of the present paper. 

Following standard definitions, a fully polynomial-time approximation scheme (FPTAS) for Z A (G) is an 
algorithm that, given as input a graph G and a parameter e > 0, outputs a number Z that satisfies 

(1 - e) • Z A (G) < Z < (1 + e) • Z A (G) (3) 

in time poly(|G|, 1/e). A fully polynomial-time randomized approximation scheme (FPRAS) is then a ran- 
domized algorithm that, with probability 1 — 5, outputs a Z satisfying ^ in time poly(|G|, 1/e, log (1/8)). 

For the Ising model, in a seminal paper [T7] Jerrum and Sinclair gave an FPRAS for Z A (G) when (3 = 7 
> 1. It was further extended to the entire region of /3j > 1 by Goldberg, Jerrum and Paterson [T3j. A two- 
state spin system is called ferromagnetic if (3j > 1 and anti-ferromagnetic if f3j < 1. The approximability 
of Z A {G) for anti- ferromagnetic systems is less well understood. Starting with counting independent sets in 
sparse graphs [ID] , the approximability of Z A (-) in bounded degree graphs is also widely studied. Significant 
progress has been made recently on the algorithmic side, and approximation algorithms for anti-ferromagnetic 
two-state spin systems have been developed in [2H [22j [191 El] > based on the technique of correlation decay 
introduced by Bandyopadhyay and Gamarnik [l] and Weitz [21]. Finally, a unified FPTAS was found [18] 
to approximate Z A (-) for all anti-ferromagnetic two-state spin systems of either bounded degree graphs or 
general graphs, when the system satisfies a uniqueness condition. 

The uniqueness condition is named for, and closely related to, phase transitions that occur for the Gibbs 
measure. It depends on not only f3 and 7 but also the degree of the underlying graph as well. Such phase 
transitions from statistical physics are believed to frequently coincide with the transitions of computational 
complexity from tractability to intractability. However, there are only very few examples where the con- 
jectured link is rigorously proved. One notable example is for the hardcore gas model (or independent set 
with f3 = and 7 = 1), for which such a conjecture was rigorously proved (for almost all degree bounds) 
both for the algorithmic side {24J and for the hardness side |23|, I12|. As discussed above |24 } 1221 IT9l I18| . for 
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general anti-ferromagnetic two-state spin systems, the algorithmic part of the conjecture has recently been 
established. In this paper, we make substantial progress on the hardness part of the conjecture. 



Our Results. For (3 and 7 : < /3, 7 < 1 except at (/?, 7) = (0, 0) or (1, 1), Goldberg, Jerrum and Paterson 
proved that the problem does not admit an FPRAS on general graphs (when there is no degree bound), 
unless NP = RP [14j . In their reduction, the degrees of the hard instances are unbounded. This is consistent 
with the uniqueness threshold conjecture. However, for any fixed /3, 7 in the unit square, the uniqueness 
condition states that there exists a finite threshold degree A(/3, 7) [22J, [19l [18] , which satisfies 

such that the system satisfies the uniqueness condition if the degree d < A(f3, r y), and the non-uniqueness 
condition if d > A(/3,7). The paper |18j gives an FPTAS for graphs with degree bounded by A(/3, 7). The 
conjectured coincidence of phase transition with hardness in complexity suggests that as soon as the degree 
of the input graph goes beyond A(/3,7), the problem becomes hard to approximate. Towards this direction 
we show that for any fixed (3, 7 in the unit square, the problem does not have an FPRAS if the degree of 
the input graph is f2(A(/3, 7)), unless NP = RP. Our hardness also holds when restricted to input graphs 
that are regular. Formally, we prove the following theorem: 

Theorem 1. There exists a positive constant h such that: Given any (3, 7 : < (3, 7 < 1 with (f3,j) 7^ (0,0), 
(1, 1) and any integer d > hj '(1 — j3~f), there is no randomized polynomial-time algorithm that approximates 
Za(G) i n d-regular graphs G with relative error e = 10~ 4 , unless NP = RP. 

Note the relation between our degree bound h/(l — /3j) and A(f3,j) from Q. 

We also make progress on (/3,7) outside the unit square. While the uniqueness condition is monotone 
inside the unit square, its behavior outside is significantly different. (See more discussions on this difference 
in Appendix [Aj) Without loss of generality, we consider the region defined by (3j < 1 with < j3 < 1 < 7. 
There is a uniqueness curve (see Figure [TJ, connecting the point (1,1) and the 7-axis. Above the curve, 
the system satisfies the uniqueness condition for any graph [19\ I18j. Hence, hardness is only possible below 
the uniqueness curve. Furthermore, when (/3, 7) is outside the unit square but below this uniqueness curve, 
there is only a finite range of degrees d for which the system does not satisfy the uniqueness condition. This 
makes it very challenging to prove a hardness result for them. Previously, the hardness was only obtained 
in p3] for a very tiny square < (3 < r] and 1 < 7 < 1 + rj where rj is roughly 10 -7 , near the point (0, 1) 
corresponding to independent set or the hardcore gas model. 

In this paper, we prove the following hardness result for (/3,7) outside the unit square: 

Theorem 2. Given (3 and 7 such that 0</?<l,7>l and /?7 < 1, let 

A' = [-l/(ln/3 + ln 7 )] and A* = [l/ln 7 ] . (5) 

When A* > 8000 A', there is no randomized polynomial-time algorithm that approximates Zj^(G) in regular 
graphs G of degree A* with relative error e = 10~ 4 , unless NP = RP. 



The new hardness region is pictured in Figure TJ^ Here the two white squares are the hardness regions 



acquired by Goldberg, Jerrum, and Paterson [14J. Beyond the uniqueness threshold, we know that FPTAS 
exists. Our hardness result, Theorem [2j applies to the region between the vertical line with 7 = 1 and the 
curve to the left of the uniqueness threshold. Let us describe the new curve in more details. Again we focus 
on the region with < (3 < 1 < 7 and f3j < 1; There is a symmetric curve when < 7 < 1 < f3. Near the 



1 The reader should be aware that, for illustration purposes, the picture is not drawn to actual scale. 
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Figure 1: The new hardness region of Theorem [2] 

point (1)1), the condition imposed by Theorem [2] is almost linear. So the new curve is roughly a line with 
slope —8000 around (1,1). When the curve approaches the line of f3 = 0, A' becomes 1 and the condition 
requires 7 to be between 1 and roughly 1 + 1/8000. 

Moreover, using a standard translation described in Appendix [B] we can generalize Theorem [T] and [2] to 
two-state spin systems with an external field. Formally, let fi > 0, we have the following two corollaries for 

Z A „(G)= £ >^<^)=°}| JJ A a{u)Mv) . 

a:V^{0,l} (u,v)£E 

Corollary 1. There exists a constant h such that, given any nonnegative /3, 7 and \x with ^7 < 1, and an 
integer d satisfying that 7 < fi d < 1 and d > izrpz > there is no randomized polynomial-time algorithm that 
approximates Zx^(G) in d-regular graphs G with relative error e = 10" 4 unless NP = RP. 

Corollary 2. Given any nonnegative (3, 7, fi and an integer d such that 

1 _1 -J— in 1 

e d < 7 • d < e d - x or e d < p ■ \i d < e d ~ x 

if d also satisfies d > 8000 [— l/(ln/3 + In 7)] , then there is no randomized polynomial-time algorithm that 
approximates Zj^^G) in d-regular graphs G with relative error e = 10 -4 unless NP = RP. 

Proof Outline. In both Theorem [T] and [2j we use the phase transition that occurs in the non-uniqueness 
region to encode a hard-to-approximate problem. This approach has been used in previous hardness proofs 
for the hardcore gas model |1U|, [2"T| I23| . To this end, we reduce the approximation of E2LIN2 to the appro- 
ximation of partition function in a two-state spin system. Here an instance of E2LIN2 consists of a set of 
variables x%, . . . , x n and a set of equations of the form x% + xj = or 1 over Z2. By |15j . it is NP-hard to 
approximate the number of satisfiable equations for E2LIN2 within any constant factor better than 11/12. 

Given an E2LIN2 instance, we use a random bipartite regular graph to encode each variable x\. Due 
to the phase transition and the fact that we are in a non-uniqueness region, each of these bipartite regular 
graphs would be in one of two types of configurations with high probability, if sampled proportional to its 
weight in the partition function. This can be used to establish a correspondence between the configurations 
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of these bipartite graphs and the assignment of the boolean variables x±, . . . ,x n . Furthermore, we also add 
external connections between the bipartite graphs according to the set of equations in the E2LIN2 instance. 
They contribute exponentially to the total weight in the partition function, according to the total number 
of equations that an assignment satisfies. Thus, a sufficiently good approximation to the partition function 
can be used to decode approximately the maximum number of equations that an assignment can satisfy. 

Our gadget is also randomly constructed. Then the probability should also be over the distribution of 
the gadgets. It is not hard to show that things work out beautifully if we simply substitute the expectation 
for the actual weight. But to make the proof rigorous, one must first obtain a sufficiently good concentration 
result. Such a result is unknown and could be very difficult to prove (assuming it is true), as it is already a 
tour-de-force in the special case for the hardcore gas model [211 E31 E2] • 

Instead we use a detour: (1) We prove a lower bound for the weights of two types of configurations we 
expect, guided by the phase transition; and (2) We prove that the total weight of other configurations is 
exponentially small compared to the lower bound with a probability exponentially close to 1. The way we 
establish the lower bound in (1) is similar to the approach by Dyer, Frieze, and Jerrum [10]. To prove (2), 
they [10] used the expectation and Markov's inequality. If we use the same approach, we could not get the 
hardness result for bounded degree graphs in the same order of the uniqueness bound. Instead, we use a 
new approach for (2). 

In fact we first prove a high concentration result for an expander property of the gadgets we use. Then 
we show that the total weight of other configurations is exponentially small, given that the gadgets satisfy 
that property. This circumvented our inability to prove a complete concentration result. But we do prove 
some limited concentration results regarding the gadget, which let us prove hardness results for degrees in 
the right order conjectured according to the uniqueness threshold. It remains open whether one can use a 
refined version of this reduction along with the proof by Sly [23] to get the exact right bound. As discussed 
in Appendix \A\ this random regular graph quite closely follows the property of phase transition in infinite 
d-ary trees, when the parameter is below or beyond the uniqueness condition. 

While the high-level idea of our proofs for both Theorem [T] and Theorem [2] are quite clear and similar, 
it remains a challenge to work out the estimation for all ranges of parameters and at the same time, make 
sure that the degree is in the same order of the uniqueness bound. To this end, technically we need to use 
very different approaches for Theorem [T] and Theorem [2j Even within Theorem [T] itself, we need to do the 
estimation differently for three different subcases. 

2 Proof of the Main Theorems 

From now on we will simply use Z{G) to denote Zx(G) or Za,u(G) whenever it is clear from the context. 

Given positive integers N and A, let H(N, A) denote the following probability distribution of A-regular 
bipartite graphs H = (U U V, E) with bipartition U, V and \U\ = \ V\ = N: here H is the union of A perfect 
matchings between U and V each selected independently and uniformly at random. (Because these perfect 
matchings are drawn independently, H may have parallel edges.) 

In both proofs of Theorem [T] and Theorem [2| we give a polynomial-time reduction from E2LIN2 to the 
approximation of Z(G). An instance of E2LIN2 consists of m equations over Z2 in n variables x±, . . . ,x n . 
Each equation has exactly two variables and is of the form X{ + Xj = b £ {0, 1}. Without loss of generality 
we may always assume m > n/2; otherwise one of the variables does not appear in any equation. Given an 
assignment S of the n variables x\, . . . ,x n , we use 6(S) to denote the number of equations that S satisfies 
and let 9* = maxs 9(S). In [15] Hastad showed that it is NP-hard to estimate 9* within any constant factor 
better than 11/12. 

Given an instance of E2LIN2 we construct a random (A + A')-regular graph G as follows, with the two 
parameters A, A' to be specified later. This construction is used in the proof of both Theorem [T] and ^1 
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Construction of G from an instance of E2LIN2. For each variable Xi, i £ [n], we let Ui and Vi 

denote two sets of dim vertices each, where di > 1 denotes the number of equations in which Xi 
appears (thus, Y^i = 2m). Moreover, Ui and Vi can be decomposed into 

Ui = U iA U • • • U U iA and V = V ijX U • • • U 

where each £/j ^ and & contains exactly to vertices. Now enumerate all the m equations in the 
E2LIN2 instance one by one. For each of the m equations do the following: 

(1) Let Xi + Xj = b £ {0, 1} denote the current equation. Assume this is the fcth time that 
Xi appears in an equation, and the Ith time that Xj appears in an equation so far, where 
k £ [di] and I £ [dj]. Denote the to vertices in U^ by {u\, . . . ,u m }, vertices in Vi t k by 
{v\, . . . , v m }, vertices in Uj^ by {u[, . . . , u' m } and vertices in Vj ; i by {v[, . . . , v' m }. All these 
vertices have degree at this moment. If b = 0, we add A' parallel edges between (u s ,v' s ) 
and (v s , u' s ), for each s £ [to]; or if b = 1, we add A' parallel edges between (u s , u' s ) and 
(v s ,v' s ), for each s £ [m]. 

By the end of this step, every vertex has degree A'. In the next step, 

(2) For each i 6 [n], we add a bipartite graph Hi = (UU Vi, Ei) drawn from H(dim, A). 

This finishes the construction, and we get a (A + A')-regular graph G with 4m 2 vertices. 

We need the following notation. Given any assignment a from V{G) to {0, 1}, we let Ui(a) denote the 
number of vertices u £ Ui with cr(n) = 0, and let Vi(a) denote the number of v £ Vi with cr{v) = 0. 

Proof of Theorem [7} Without loss of generality, we assume < (3 < 7 < 1 . We can also assume that f3 > 
since the tight hardness to the exact uniqueness bound for (3 = has been shown in [19] by generalizing the 
tight hardness result for the hardcore model |23l H~2] . 

Given any assignment S of the n variables, we let Z(G, S) denote the sum of lj(G, a) over assignments 
a : V(G) — > {0, 1} that satisfy for each i £ [n], 

Ui(a) < V(a) if Xi = in S; or Ufa) > Vfa) if x t = 1 in S. (6) 

By definition, we have Z(G, S) < Z(G) < J2s %(G, S). We prove the following key lemma in Section [3] 

Lemma 1. There exists a positive constant h with the following property: for any f3 and 7: 0</3<7<l 
with (/3,7) 7^ (1, 1) and for any A* > /i/(l — (3j), there are D > 1, C > and positive integers A and A' 
with A + A' = A* ; that satisfy the following property: given any input instance of E2LIN2 with n variables 
xi, . . . ,x n and m equations, except for probability < exp(— f2(m)), the A* -regular graph G constructed with 
parameters A and A' satisfies 

C vn? . D m0(S) < z{ q^ S j < C m 2 . (e(S)+0.03m) ^ j Qr any assignment S Q f the n varia bles. (7) 

Given /3,7 and A*, we let C, D, A and A' denote the constants that satisfy the condition in Lemma [TJ 
then given an input instance of E2LIN2, ([7]) holds with probability 1 — exp(— Q(m)). 

Now assume holds. We use 9* to denote the maximum number of consistent equations and use S* to 
denote an assignment that satisfies 9* equations. We also use Y to denote an estimate of Z = Z(G), where 
\Y/Z - 1| < e = 1(T 4 . From Q and Z(G, S*) < Z{G) < J2s Z ( G > S )> we § et 

(1 + e) • 2 n • C m2 ■ D m ^* +omm ^ > (1 + e) • Z > Y > (1 - e) • Z > (1 - e) • C" 1 ' • D md * (8) 
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Using Y, we set 

, _ InY - ln(l + e) - nln2 - m 2 ln(7 - 0.03m 2 InD 

mlnD 

and we get Y' < 9* since InD > 0. We finish the proof by showing that Y' > (11/12) • 9*. By Q we get 

, > ln(l + e) - ln(l - e) + nln2 + 0.03m 2 InD 

mliiD 

As 9* > m/2 and m > n/2, when m is large enough, Y' > (11/12) • 9* and the theorem is proven. □ 

Next, we prove Theorem [2] 

Proof of Theorem^ For /3,7 with < f3 < 1 < 7 and /37 < 1, let A' and A* be the two positive integers 
defined in (jsj) which satisfy A* > 8000 A'. We set A = A* - A'. Given any input instance of E2LIN2 with 
n variables and m equations, we use G to denote the A*-regular graph constructed using A and A'. 

First of all, we show that to get a good approximation of Z(G), with high probability it suffices to sum 
u(G,a) only over assignments o~ that satisfy the following condition: 



min (Ui(a), Vi{a) j < Xckm, for all i G [n], where A = 9 x 10~ 5 . (9) 

We use X to denote the set of such assignments. Formally, we prove the following key lemma in Section |4j 

Lemma 2. Let G be the graph constructed from an instance of E2LIN2 with n variables x\,X2, ... ,x n and 
m equations, with parameters A and A'. Then with probability 1 — exp(— ^(m 1 / 3 )), it satisfies 

£>(G» < Z(G) < (l + o(l)) • J>(G». (10) 
o-es o-es 

Next, given an assignment S over the n variables we use Zy,{G, S) to denote the sum of co(G, a) over all 
assignments a 6 X that satisfy ^ for all i G [n]. We prove the following lemma: 

Lemma 3. There are C > and D > 1 satisfying the following property: given an instance of E2LIN2 with 
n variables and m equations, the A* -regular graph G constructed with parameters A and A' satisfies 

C m2 ■ D md ^ < Ze(G, S) < C m2 ■ D m ( e ^ +0Mm ) , for any assignment S of the n variables. (11) 

Proof. Let S be an assignment over the n variables, and we use the same lower bound 

Z S (G, S) > Z*(G, S) = C rn2 ■ D me ^ 



with Z*(G,S) defined in (15), C and D defined in (16). It is a lower bound for Z-£(G,S) because every a 
that satisfies ( |14[ ) is in X by definition. Since j3j < 1, we have C > and D > 1. 

Now we give an upper bound for Z-s(G, S). For each a in the sum Z*(G, S), we use Q a to denote the 
following set of assignments a 1 in the sum of Zy,{G, S): For each i £ [n], 

1. If Xi = in S, then <t' agrees with a on and Ui{a') < min { Vi(cr), Ad^m}, while Ui(a) = 0; 

2. If Xi = 1 in 5", then a' agrees with a on U\ and V^c') < min { Ui(a), A^m}, while V^(cr) = 0. 

It is easy to show that {Qa} is a partition of the assignments in Z^(G, S). Moreover, as in ( |21[ ) we have: 

\QJ <m 2n -e 2H ^ M 
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Every a' in Q a has weight u(G,a') < ui(G, a) because nipping a bit from 1 to cannot increase its weight. 
Finally, we get the following bound for Z^(G, S): 

Zx(G, S) < Z*{G, S) ■ m 2n ■ e 2H{x)M 

To finish the proof, we plug in A = 9 x 1CT 5 and H(X) < 0.000929 to compare 2H(X) with 0.04 In D. Also 



recall the definition of D in (16). From the assumption that 

(/3 7 ) A ' < 1/e and 7 A+A ' < cy < l.OOle 
we get lnL> > 0.04673 and 0.04 In L> > 2H(X). (Here 7 < 1.001 because by the condition of Theorem [Tj 

e > 7 A * _1 and A* > 8000 A' > 8000 

and thus, 7 < 1.001.) This finishes the proof of the lemma. □ 

Let 8* > m/2 denote the maximum number of consistent equations, and let S* denote an assignment 
that satisfies 9* equations. From these two lemmas, we have with high probability that 

C™ 2 ■ D md * < Z^(G, S*) < Z(G) < (1 + o(l)) • 52 s Z s(G, S) < (l + o(l)) • 2 n • C™ 2 ■ D m ( e *+ 0Mm ) 

Theorem [2] then follows from the same argument used in the proof of Theorem [TJ □ 

3 Proof of Lemma [T] 

Recall that /3 and 7 satisfy < (3 < 7 < 1 and (/3, 7) ^ (1, 1). 

Given any instance of E2LIN2 with n variables x%, . . . ,x n and m equations, the (A + A')-regular graph 
G we construct consists of A-regular bipartite graphs Hi, i £ [n], and edges between them. For each Hi in 
G, we use Ui U V% to denote its vertex set with \Ui\ = \Vi\ = dim. Recall that e = 10~ 4 . We say A C Ui is 
big if |^4| > e|f/j|, and B C Vi is big if \B\ > e\Vi\. We also use E(Hi, A, B) to denote the number of edges 
between A and B in Hi. Let K = 48/e 2 = 48 x 10 8 . We prove the following lemma: 

Lemma 4. When A > K , with probability 1 — exp(— £l(m)), the graph G we construct satisfies 

E(Hi, A, B) > ^ilL^i ) f or a ll % e [n] and for all big A C Ui and big B C V;. (12) 
4cijm 

Proof. We let re denote dim and let d denote A. Let U, V be two sets of re vertices. Let H be a random d- 
regular bipartite graph generated by independently picking d perfect matchings between U and V uniformly 
at random. We prove the following lemma. Lemma [4] then follows from Lemma [5] using the union bound. 

Lemma 5. Given any AC.U and BCF with \A\ = an and \B\ = bn, where b > a > 10 -4 , we have 



Pr 



the number of edges between A and B in H < abdn/A < 2 cdn , where c = 1/ (16 • 10 8 ). 



Proof. The intuition is that H is drawn from a distribution that is very close to G(n, d). 

We use ui, . . . , u an to denote the vertices in A. For each k £ [d] and i E [an], we use Yj^ to denote the 
random {0, l}-variable such that Yk t i = 1 if U{ is matched with a vertex in B in the /cth perfect matching; 
and is otherwise. From this, we have 
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Pr 



the number of edges between A and B in H < abdn/A 



Pr 



Yfc^ < abdn/4 

k£[d], i£ [an] 



However, the variables Yfc,!, . . . , Yk^n are clearly not independent. 

To deal with this issue, we introduce the following independent random {0, l}-variables X k ^ for every 
k G [d] and i G [an]: X k ^ = 1 with probability pi = (bn — i + l)/n and = with probability 1 — /jj. It 
is easy to see that Y k; i dominates X kt f For all k and z and for all y k i, ■ ■ ■ , Vk,i-i G {0, 1}, we have 



Pr 



Yhl = Vk.i, ■ ■ ■ ,Yk,i-l = Vk,i-l 



bn - (y kjl H h y k , 



i-lj 



n 



> Pi = Pr 



Next we define, for each k, a sequence of {0, l}-random variables Z k ^. They are highly correlated with the 
Xk/s but have the same joint distribution as {Y k> i : i G [an]}: 

1. First, Z k ,i = X k>1 . Thus, Pr [Z k .x = l] = Pl = b = Pr [Y kA = l] . 

2. For Zk : 2- If X kj 2 = 1) then Zk t 2 = 1> If ^k,2 = (which happens with probability 1 — P2), then we 
set Zk,2 randomly as follows: Suppose Z k> \ = z kt \ G {0, 1}, then set Z ki 2 = 1 with probability 



(bn - z kl )/n - p 2 



1 - P2 

and set Zk 2 = with probability 1 — w. Note that we have 



Pr 



Z 



k.2 



1 



ZkA — Z k A 



p 2 + w(l - P2) = {bn - z k ,i)/n. 



Thus, \Z ki i, Zk^} has the same joint distribution as lfc,2}- We also have Z kj 2 > ^k,2- 

3. For Zk,i from i = 3 to an: If -X^j = 1, then we set Z k j = 1; If X^i = (which happens with 
probability 1 — pi), then we set Z k ^ randomly as follows: Suppose Zk j = z k j G {0, 1} for each 
j G [i — 1], then set Z k i = 1 with probability 



w 



1 (bn-(z k>1 -\ \- Zk,i-i))/n- pi 



1 - Pi 



and set Z k i = with probability 1 — w' . Then note that 



Pr 



Z, 



k.i 



Zk,i = Zk,i, ■■■ , Zk,i-i = Zk,i-i = Pi + w'(l - pi) = (bn - (z k> i H V z k)i -i)) j n 



Thus, {Z kt i, ■ ■ ■ , Z k ,i) has the same joint distribution as • • • , We also have Z k ^ > X k ^. 

From the construction above it is clear that the joint distribution of {Zk,i : i G [on]} is the same as that of 
{Yk,i '■ i G [an]}. Because Z k> i > Xk^ for all i G [an], we have 



Pr [ Y k,i > abdn/A = Pr [ ^ Z k ,i > abdn/A > Pr [ ^ X kji > abdn/4 



(13) 



By (13), it now suffices to prove an upper bound for 
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We can now use the Chernoff bound. The expectation is \i = i ELX"^] = d • ^ pi > abdn/2. 
By the Chernoff bound (and setting 5 = 1/2), we have 

Pr[ Y, Xk ' i - abdn / A <Pr[5^X fc)i < (l-5)n <exp(-^ 2 / 2 )- 

The lemma then follows from 

, r2/ , / abdn / 1\ 2 l\ ( 10~ 8 dn/l\ 2 l\ 
exp(-^/2)<exp^-— • - j < exp ^ 

by setting c to be l/(16 • 10 8 ) . □ 
Using Lemma [5| we can then choose a large enough d and apply the union bound on A and B. □ 

To finish the proof of Lemma [TJ we divide (/3,7) into three cases. For each case we show there exists a 
large enough constant h with the following property: for all (/3, 7) of this case, and for all A* > h/(l — /S7), 
there are C > 0, D > 1 and positive integers A > K and A' > 1 with A + A' = A*, such that holds 
whenever G satisfies (|12l). The lemma then follows by taking the maximum of the three /i's. 



Let L = 12/e = K/4 = 12 x 10 . In the rest of the proof, we assume G satisfies (12) and let M = m . 



3.1 Case 1: < (5 < 1/2 and /3 < 7 L 

We set h to be a large enough constant so that h/(l — fij) > 7(L + 1). Given a A* > h/(l — ^7), we then 
set A = [LA*/(L + 1)J > K and A' = \A*/(L + 1)] > 7, where A + A' = A* and LA' > A > L(A' - 1). 

Let S be an assignment over the n variables x±, . . . ,x n , then we start with a lower bound Z*(G, S) for 
Z(G, S). To this end, we consider the sum of oj{G, a) over all assignments a that satisfy for each i £ [n]: 

Ui(a) = if Xi = in 5; and V^(cr) = otherwise. (14) 

Denote this sum by Z*(G, S). It is clearly a lower bound for Z(G, S), and is exactly equal to 

Z*(G, S) = (l + 2 7 A + A ' + ^^') me(S) • (V' 7 A' + 27 A + A' + 7 2A + 2A^(™-^)) 



md(S) 

" - " - 2 ~ ^ + - ^ ) ■ 1 ^>"; 27A+A / +72A+2A , 1 • (15) 



\ M I 1 1 o ,A+A' 1 ..2A+2A' 
jA'.A' , o-A+A' , 2A+2A'V W I I + 27 +7 



Setting C and D appropriately, we have Z*(G, S) = C M ■ D m8 ^ , where 



1 , 2 ^ A + A ' 4- ^,2A+2A' 
2A+2A' ^ n„„j n ' 1 ^ 7 



C = 0* 7 A + 2 7 ™ + 7 ^ + ^ > and /J = j ^ A+A , + yA+2A/ > 1 (16) 

since (/3, 7) 7^ (0,0), (1, 1). It is also easy to give a lower bound of 8/7 for 7J because the difference between 
the numerator and the denominator is 1 — /3 A 7^ > 1/2 as /3 < 1/2; and the denominator of D is < 7/2. 
Next, to give an upper bound for Z(G, S), we consider the sum of oj(G, a) over a that satisfies ^ and 

Ui(a) < edivn if xi = 0; and Vi(a) < edim if = 1 (17) 

for every i 6 [n]. We show that this sum is indeed a good approximation of Z(G, S): 
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Z(G,S)<(l + o(l)) Yl ^ G ' CT ) 

cr satisfies (JsJ and ( |17| for all i 



(18) 



To show (18) we randomly draw an assignment a from those appear in the sum Z(G, S) with probability 
proportional to co(G, a), and it suffices to show that the probability that a satisfies (17) for all i is 1 — o(l). 
This then follows from the following lemma and the union bound: 



Lemma 6. For any i 6 [n], the probability that a violates (17) is at most exp(-dim). 

Proof. Without loss of generality, we assume Xi = in S. 

Pick any partial assignment a' over all vertices of G except those of Hi which satisfies ^ for all j 7^ i. 
To prove the lemma, it suffices to show that the sum of co(G, a) over all assignments a that are consistent 
with a' and satisfy ^ for i but violate (17) at i is exponentially smaller than uj(G, a*), where we use a* to 
denote the unique assignment that is consistent with a' and satisfies Ui(cr*) = and Vi{a*) = dim. 

To this end, we let u(a') denote the product of the edge weights in a' over all edges in G except those 
have at least one vertex in Hi. Then it is easy to give the following lower bound for uj(G, a*): 



A'im 



(19) 



On the other hand, for any assignment a that satisfies ([6]) (hence, a satisfies Vi(a) > Ui(a), as we assumed 
Xi = 0), is consistent with a' but violates (17), we must have 



u(G,a) < u(a') • 



<E 2 Adim/4 



(20) 



It follows from the assumption of (12) and Vi{a) > Ui(a) > ediin as Xi = and a violates (17). 
Plugging in A and A', we have 



cu(G,a*) 07(0 • (07) 



A' dim 



> 







2A' 



dim 



u(G,a) ~ w(cr') ./^ 2 a<W4 - I 03(A'-i) 



1 



A'-3 



> 2 



3dim 



The lemma follows because the number of a that is consistent with a' but violates (17) is at most 2 2d ' m . □ 



We continue with (18). For each a that satisfies (14), let T a denote the following set of assignments a': 
(i) a' satisfies (17) for all i\ and (ii) for each i, a' agrees with a over Vi if Xi = 0; and agrees with a over Ui 
if Xi = 1. It is clear that {T a } is a partition of the assignments that satisfy (17) for all i £ [n]. It is easy to 
check that for any a that satisfies (14), we know exactly the cardinality of \T a \: 

[edim\ 





<n((N,.»j + i)( L *: 



ien 



< m 



2/i 



TJ e H{e)d i m = m 2n. e 2H{e)M ^ 



IS n 



where -ff(e) ~ 0.00102. For any a 1 £ T CT , we also have 



UJ 



(G,a') <u>(G,a)- (l/ 7 



A+A' 



2eM 



(22) 



because for any assignment, switching the value of a vertex from 1 to can improve co(G, a') by at most a 
factor of 1/7 A+A '. Finally, by combining ( 18 ), (21 ) and (22 ) we get the following upper bound for Z(G, S): 



Z(G, S) < (1 + o(l)) • m 2n ■ e 2H(e)M ■ Z* (G, S) ■ (l/ 7 



A+A' 



2eM 



(23) 
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To finish the proof and show Q: Z(G, S) < Z*(G, S) • D 03M , we also need to compare D with 1/7 A+A '. 
Since /3 < 7 L , we have /3 A ' < j LA ' < 7 A . It follows from the definition of D that D > 1/(47 A+A '). Then 
(|7j) follows directly from (23) by plugging in D > 8/7 and e = 1CP 4 . 



3.2 Case 2: (3 > 1/2 and /3 < 7 L 

We set /i, A and A' as follows. We pick h to be a large enough constant so that for any A* > h/{l — /S7), 
A = L^A*/(L + 1)J and A' = \A*/(L + 1)J satisfy A' > 7 and (f3j) A> < 2~ 12 . Since 7 > 0, we also have 
/3 A ' < 1/64. By the definition of A and A', we have LA' > A > L(A' - 1). 

The proof follows the same flow as that for Case 1. First of all, we use the same lower bound in (15): 

Z(G, S) > Z*(G, S) = C M ■ D m9 ^ 



where Z*(G, S) is the sum of co(G,a) over a that satisfies (14). C > and D > 1 are defined as in (16). It 
then follows from (/?7) A ' < 2~ 12 that D > 4/(3 + 2" 12 ) » 4/3. 



We also use the same upper bound argument (18) for Z(G,S): 

Z(G,S)< (l + o(l)) J] w(G,<7). 

(T that satisfies | |17| l 



(24) 



To prove the same statement as in Lemma [6j pick an i G [n] and any partial assignment a' over all vertices 
of G except those of Hi, satisfying ^ for all j / i. Without loss of generality assume Xi = in S. We use 
a* to denote the assignment that is consistent with a' and satisfies Ui(a*) = and Vi(a*) = dim. We also 
use a to denote any assignment that is consistent with a', satisfies ^ in Hi, but violates (17). Then from 
(19) and (20), we have 



A' dim 



m(G,a*) > (fry) 

W(G» ~ /^ 2 Ad»m/4 




dim 



> 



1 

/3A'/2 



dim 



> 2 



3dim 



The upper bound ( 24 ) then follows directly. To finish the proof we define T a similarly for each assignment 



a that satisfies (17). Combining (j24j), (21) and (22), we obtain the same upper bound p3[ ) for Z(G,S). It 
also follows from (3 < j L that D > 1/(47 A+A '). Then Q is proven by plugging in D ~ 4/3 and e = 10~ 4 . 

3.3 Case 3: l3 > 7 L 

For this case, we need to use a different estimation for Z(G, S). 

We start by setting A' and A. Let h be a large enough constant such that for any A* > h/(l — (3j), 



A 



L(L + 1) ■ A* 
L(L + 1) + 1 



and A' 



A* 



L(L + 1) + 1 



satisfy A' > 1 and (/?7) A ' < 1/4. It follows from the definition that A* = A + A' and A > L(L + 1)A'. 

Given an assignment S over the n variables x\, . . . ,x n , we use a to denote the unique assignment with 
Ui(a) = and Vi(&) = dim when X{ = 0; Ui(a) = dim and Vi(a) = when x\ = 1, for all i G [n]. Then 



Z(G,S)>u{G,a) = (0y) 



A'm(m-6(S)) 



A' 



M 



m8{S) 



(25) 



Setting C = (/3 7 ) A ' and D = 1/C, we get Z(G, 5) > C M ■ D md ^ s \ with D > 4 and G > 0. 



11 



Next, to give an upper bound for Z(G, S), we consider the sum of u)(G, a) over a that satisfies 



Ui{cr) < e ■ dim and Vi{<j) > (1 — e) ■ dim, when Xi = in S; 
Ui(cr) > (1 — e) • dim and Vi(a) < e ■ dim, when Xi = 1 in 5". 

for every i E [n]. We show that this sum is a good approximation of Z(G, S): 

Z(G,S)< (l + o(l)) £ u;(G» 

cr that satisfies J26l and J27I 



(26) 
(27) 



(28) 



To prove (28), we randomly draw a a from those appear in the sum Z(G, S) with probability proportional 
to uj{G,a), and show that the probability that a violates (26) or (27) is exponentially small. 

For this purpose we prove the same statement as in Lemma 6. Pick any i £ [n] and assume Xi = in 
S without loss of generality (the proof below handles (26), and (27) is by the same argument.) Let a' be 
any partial assignment over vertices of G except those of Hi, and satisfying ^ for all j ^ i. We use a* to 
denote the unique assignment that is consistent with a' and satisfies Ui(cr*) = and Vi{a*) = dim. We let 
a denote any assignment that is consistent with a' but violates (26) in Hi. Then we get 

dim 



u(G,a*) Qj(a') • (/3 7 ) AV ' m ( ^ L+ ^ A ' 



3(L+1)A' 



7 



1 



2(L+1)A' 



dim 



> 2 



Adim 



Here the first inequality follows from (12) and the fact that, since a violates (26), either 

Vi(a) > Ui(a) > e ■ dim or Ui(a) < Vi(a) < (1 - e) • dim 



(29) 



and the last inequality follows from 1/4 > (/?7) A ' > r y^ L+1 ^ A ' . This proves (28). 



Moreover, the number of a that satisfies (26) and (27) for all i 6 [n] can be bounded by m 4n • e iH ( e ) M ; 



in the same way as (21). And for any a that satisfies both (26) and (27), we also have 

,2eA'M 



u(G,a)<u(G,a)/(p 1 y 



(30) 



This is because, to obtain a from a, each time we flip a vertex from to 1, the weight increases by a factor 
of at most 1/7^'; each time we flip a vertex from 1 to 0, the weight increases by a factor of at most 1//3 A '. 
Finally, combining (28), (29) and ( pO] ), we get 

Z(G, S) < (1 + o(l)) • m 4n • e m ^ M ■ u(G, a) ■ D 2eM 

Then aln follows immediately by plugging in D > 2 and e = 10~ 4 . This finishes the proof of the lemma. 



4 Proof of Lemma [2] 

Recall that f3 and 7 satisfy /3,7:0</3<l<7 and /?7 < 1. Let A' and A* be the two positive integers 
defined in Theorem [2J with A* > 8000A'. From their definitions, we have (/3 7 ) A ' < 1/e and 7 A * > e. Set 
A = A* - A' > 7999A' > 7999. By the definition of A*, we have e > 7 A *^ > 7 A and thus, 7 < 1.001. 

Given an instance of E2LIN2 with n variables x\, . . . ,x n and m equations, we use G to denote the A*- 
regular graph constructed with parameters A and A', where A* = A + A'. We let Hi denote the bipartite 
graph in G that corresponds to Xi and use Ui U Vi to denote its vertices, with \Ui\ = \ Vi\ = d^m. 

Before working on G and Hi, we start by proving a property that a bipartite graph sampled from the 
distribution H(N,A) satisfies with high probability. Let H be a bipartite graph drawn from T~L(N,A) for 
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some N > 1 and A denned above, with 2N vertices U UV. We also use p : U U V — > {0, 1} to denote an 
assignment and call it an (a, 6)-assignment for some a, b £ T/v, where Tv = {0, 1/N, 2/N, . . . , (N — 1)/N, 1} 
if \ueU : p(u) = 1 = aN and | v G V : p(v) = 1 = bN . We also use X/v(a, b), where a, 6 G T/v, to denote 
the set of all such (a, 6)-assignments, and let 



P€In (a,b) 

with A' defined above. We are interested in the expectation of Z a ^{H) when min(a, b) > A = 9 x 1CT 5 : 
Lemma 7. For large enough N and a,b G T/v suc/i that min(a, b) > X, we have 



(31) 



E 



Z a ,b(H) <exp(l.21-7V). 



Proof. We recall the definition 



Z atb (H)= E ^(^p)-7 A ' (2 - a - b)iV . 

pel N (a,b) 



We want to compute the expectation 



E = E 



Z a ,b{H) 



Because the distribution H(N,A) is symmetric for any p £ I^(a,b), each term in the summation Z a ^(H) 
has exactly the same expectation. Thus, we have 



E = ^'(2-a-b)N . 



N\ ( N 



aN \bN 



E 



QkN -,(l-a-b+k)N (bN\ /(l-b)N\ 
P ' \kN) \(a-h)N) 



, fceTjv, a+fe— l<fc<min(a,6 



Since we only care about the exponent of E and the summation of k is only over a linear number of terms we 
can replace the summation by maximum without changing the leading term of the exponent. Therefore we 
conclude that the coefficient of iV in the exponent of the expectation E = exp(\I/(a, b)N) is of the following 
form, in which the function H(x) = — xlnx — (1 — x) ln(l — x): Let c = 8000, we have 



*(a,6) 



max 

k 



max 

k 



< max 

k 



< max 

k 



(2- a- 6) A' In 7 + #(a) + H(b) + A^fc ln/3 + (1 - a - 6 + k) In 7 + 6i7(|) + (1 - fc)^!^) - # (a)) 
A' In 7 + (1 - a - 6)(A + A) In 7 + fl-(a) + if(6) + A^ln(/3 7 ) + bH(*) + (1 - - tf(a)) 

^ + (1 - a - b) + (a) + if(6) + (c - 1) (fc A' ln(/3 7 ) + + (1 - (f^f) " # (a)) 

^ I + (l-a-6)^ I +^(a) + ff(6) + ( C -l)(-fc + 6ff(|) + (l-6)^(^|)-^(a)) 



For the last formula, we can use Mathematica to verify that its value is < 1.21 when min(a, b) > 9 x 10 5 . □ 



13 



For now, by Lemma [7] we can impose the following condition on the graph G constructed from the input 
instance of E2LIN2: 



For all i G [n] and all a, b G T^ im with min(a, b) > X, Z a ^(Hi) < exp(l.22 • dim). 



(32) 



Using Lemma [7j Markov's inequality and the union bound, it is easy to show that G satisfies this condition 
with probability 1 — exp(— Q(m)). In the rest of the proof we show that G satisfies (10) whenever it satisfies 
(32). Lemma [2] then follows immediately. 



We assume that G satisfies (32). To prove (10), we randomly sample an assignment a with probability 



proportional to u(G,a). (10) follows if we can prove that a satisfies ^ with probability 1 — o(l). For this 
purpose we need the following lemmas that give us properties that a satisfies with high probability. Given 
any set L of vertices in G, we let N a (L) = {v G L : o~(v) = 0}. Also recall the definition of Uik and Vi^ in 
the construction of G. Let X{ be a variable and let k G [di] , then 

Lemma 8. Let a be an assignment drawn according to its weight. For any i G [n] and any k G [di], except 
for probability exp(— ^(m 1 / 3 )), we have 



\N a (Ui, 



< 



1 + e 



1 + m- 1 ^) ■ \ U, 



Proof. Pick any partial assignment a' over vertices of G except those in Un~- Conditioned on a' , it is easy 
to see that the values of vertices in are independent. Each vertex in Ui k has A + A' neighbors, each 
of which contributes a vertex weight of either /3 or 1 if it is assigned 0, and either 1 or 7 if it is assigned 1. 
Since 7 < 1//3, the total weight for assignment 1 is at least 7 A + A ' > e times the weight for assignment 0. 
The lemma follows from the Chernoff bound. □ 

Given an assignment a, we use cr, to denote its restriction over vertices in Hi and <7_j to denote its partial 
assignment over vertices in G except Hi. We let M a _ i (Ui) denote the subset of Ui whose unique neighbor 
outside of Hi is assigned 1. Using Lemma [8] and the union bound, we have 

Corollary 3. Let a be an assignment drawn according to its weight. Except for probability exp(— Sl(m 1//3 )) 



1^(^)1 > 



1 + e 



0(m" 1/3 ) ) • \Ui\, for alii G [n]. 



(33) 



It is also clear that Lemma M and Corollary [3] also hold for and Vi, respectively, by symmetry. Now 
we are ready to prove Lemma Let a = (o"j,cj_j) be an assignment drawn from this distribution. Recall 
the definition of X below ^ . Then by Corollary [3] we have 

Pr[cr $ S] < exp(-ft(m 1/3 )) + Pr a S cr_i satisfies ([33]) for both Ui and Vi and for all i G [n] (34) 



To prove an upper bound for (34) we fix cr_j to be any partial assignment over the vertices of G except 



those of Hi, which satisfies (33) for both Ui and Vi. Then it suffices to prove that the sum of cj(G, a) over 
all g G S that are consistent with <r_j, denoted by Z\, is exponentially larger than the sum of u>(G, a) over 
all a £ X that are consistent with cr_j, denoted by Z<i- 

Let u(o~-i) denote the product of the edge weights in o~-i over all edges that have no vertex in Hi. By 



the definition of Z a ^(H) in (31), we have 
Z 2 < w(cr_i) 



^ Z a:b (Hi) < uj(a-i) ■ (dim) 2 ■ exp(l.22 • dim), 

a,beT d , m :a,6>A 



(35) 
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where the second inequality follows from (|32j). To prove a lower bound for Z\, we let 

L = | M ff _ t (Ui)\ and R = \ M a _ t (Vi) \ . 
Consider all assignments a that are consistent with <r_j and Ui(a) = 0. This gives us 

Z X > w(a_ f ) • 7 A ' L • (1 + 7 A+A ') R • W A ' + 7 A f> m - R . 



By plugging in > e, 7 A > e 7999 / 8000 ) a g we ll as the lower bound valid for R in (33), we get 

Zi > w(<r_ f ) • exp(l. 22897 • ckm), 



and the lemma follows from ( 35 ) . 
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A Uniqueness 

Fix /3,7 > such that f3j < 1. We define the uniqueness condition and discuss some of its properties. 
Definition 1. let x be the positive fixed point of the following function 



This condition specifies whether the Gibbs measure of the system on the (d + l)-infinite regular tree is 
unique or not. 

For a fixed tuple of parameters (/3, 7,/i), if (/?, 7) lies inside of the unit square, i.e. /3, 7 : < (3, 7 < 1, 
as d increases, the uniqueness condition will eventually fail. In this region, the monotonicity property with 
respect to d holds. That is, there exists some threshold A depending on ((3, 7, /i) such that, for any d > A, 
the uniqueness condition does not hold. 




We say that a tuple ( j3, 7, fi) exhibits uniqueness on d-regular graphs if 



(x + 7) d+1 



(0x + l)(x + 7) 



d(l - pj)x 



< 1. 
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When one of /3,7 is larger than 1, to satisfy the uniqueness condition, the dependence between 7, u) 
and d gets tricky. There exists a boundary such that, beyond it the uniqueness condition always holds. On 
the other hand, within the boundary, when d increases, the uniqueness condition will eventually fail. But if 
d gets even larger, the uniqueness condition will hold again. 

A.l Random Regular Bipartite Graphs 

A classical gadget to use is random regular bipartite graphs by combining d perfect matchings. This is also 
the gadget we used in this paper. One intuitive reason why this gadget is good is because of its tree-like 
local structure. More formally, we can show that the expected behavior of this gadget undergoes a phase 
transition when the parameters of the system go across the uniqueness boundary. Let Z a ^ be the expected 
weight summing over only subsets of size an and bn assigning on each side respectively. If the system 
satisfies the uniqueness condition, then the maximum of Z a ^ is achieved at the single point (p*,p*)', while 
if the system does not satisfy the uniqueness condition, Z a ^ achieves its maximum at two points (p + ,p~), 
(p~~ ,p + ), where p + > p~. So the idea is to use these two maximum points to encode two states (a {0, 1} 
assignment of a variable, two parts of a cut, and so on), and reduce another problem to this. 

The difficulty here is that the connection described above is only proved with respect to expectations. 
To make the reduction go through, we need to show that it holds with high probability, when we randomly 
choose a fixed gadget. Proving such a high concentration result, however, is not easy. For the special case 
of the hardcore model (/3 = 0, 7 = 1) such a high concentration result was almost obtained after a sequence 
of work by a careful analysis of its second moment |21| 1231 [T2] . Proving such a result for general two-state 
spin systems seems very difficult. 

Instead we observed that as d increases, p + approaches 1 and p~ approaches 0. Hence, as long as d is 
large enough, we still manage to get some exponential gap between the weights of points near (0, 1) or (1, 0) 
and everywhere else that holds with high probability. We further argue that to achieve this gap, the degree 
d differs with the desired uniqueness threshold by only a constant factor. 

A. 2 The order of the threshold degree 

Given the definition of the uniqueness condition, we are not able to give a closed form for the boundary d 
of uniqueness and non-uniqueness in terms of f3, 7 and u in general. But some properties were known and 
we summarize them as follows. The following lemma is from [22\ and [18J: 

Lemma 9. Assume /S, 7 > and Pj < 1. ///3,7 and d satisfy 




or equivalently, 




the system always satisfies the uniqueness condition for any external field u. 



Because 




l-/3 7 ' 



1 
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For d > (1 + V?7)/(l - \//37), we define 

-1-07 + d(l - Pi) - ^(-1 " Jl + d(l - Pi)) 2 ~ 4/37 



x 2 (d) 



2/3 

-1 - /3 7 + d{\ - pi) + ^{-l- Pi + d{l- Pi)Y-APi 
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which are the two positive roots of the following equation 

d(l - Pi)x 
(Px + l)(x + i) ~ 

The following lemma is from [18] . 



Lemma 10. Ifj > P > 0, pi < 1 and \J~Pi < [d— l)/(d+l), then the system described by (P,j,fi) exhibits 
uniqueness on d-regular graphs if and only if < ^i(d) or [i > /-12(d), where 

Hi(d) = Xi(d) ( Xl f + 7 ) , fori = land2. 
\Pxi{d) + 1/ 

When < P, 7 < 1 are treated as constants and [i approaches zero or infinite, we get the dependence of 
d in terms of /i as follows: 

• When d is large, 21(d) is very small, n±(d) is in the order of d 7 d . To get non-uniqueness 
for very small /i, we need that /i > /ii(d). This gives the bound log ///log 7. 

• When d is large, X2(d) is very large, /12(d) is in the order of d/p d . To get non-uniqueness 
for very large /x, we need that fi < [12(d). This gives the bound — log///log/3. 

Therefore, the dependence of d in terms of \x in Corollary [T] is also tight up to a constant factor given 
that < P, 7 < 1 are treated as constants. 



B Spin Systems with External Field 

It is easy to verify that Zjs L ^(G) can be written as 

o-:V-S>{0,1} 

where we use s(a) to denote the number of v € V with o~(v) = 0; to(o~) to denote the number of (u, v) G -E 
with <r(w) = = 0; and t\(a) to denote the number of (u,v) £ E with a(u) = a(v) = 1. 

Let ^(c) denote the number of edges whose two ends are assigned different spin states in a. For a regular 
graph of degree d, we have d ■ s(a) = 2to(a) + £2(0"), and to(a) + h(a) + £2(0") = Thus we can get 

2to(fj) + t 2 (^) |#|+t (<7)-*i(a) 
S(CJ) = d = d 

and rewrite Za_^(G) as 

ZaA g)=^ y: , tji ^-P^-^=^ e (^y° ia) ■ U) tliCT) 

a:V^{0,l} a:V^{0,l} J 
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The global factor \i a can be computed in polynomial time and the summation part can be considered as 
the partition function on the same graph with the following new parameters 

(/3',7V)= f^'^' 1 ) 

This is a two-state spin system without external field. We can apply Theorem [T] directly to get Corollary [TJ 
To apply Theorem [2] we also need to further impose that the degree satisfies a certain relationship with the 
weight after transformation. This explains the conditions specified in Corollary [2| 
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